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Abstract 

Relativistic mean-field calculations of multi-A" hypernuclei are performed by adding K~ mesons 
to particle-stable configurations of nucleons, A and 3 hyperons. For a given hypernuclear core, the 
calculated K separation energy saturates with the number of K mesons for more than roughly 
10 mesons, with bounded from above by 200 MeV. The associated baryonic densities saturate at 
values 2 — 3 times nuclear-matter density within a small region where the ^T-meson densities peak, 
similarly to what was found for multi-^ nuclei. The calculations demonstrate that particle-stable 
multistrange {N, A, 3} configurations are stable against strong-interaction conversions A — > NK 
and 3 — > NKK, confirming and strengthening the conclusion that kaon condensation is unlikely 
to occur in strong-interaction self-bound strange hadronic matter. 
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I. INTRODUCTION 



Quasibound nuclear states of K mesons have been studied by us recently in a series of 



articles 



i l2l, li| , using a self-consistent extension of nuclear relativistic mean- field (RMF) 
models. References [1, 

as 

focused on the widths expected for K quasibound states, 
particularly in the range of K separation energy ~ 100 — 150 MeV deemed relevant 
from K~-aXora phenomenology f^llj] and from the KEK-PS E548 12 C(K~ , N) missing-mass 
spectra [6|] that suggest values of Re V^(po) ~ —(150 — 200) MeV. Such deep potentials are 
not reproduced at present by chirally based approaches that yield values of Re Vg-(po) °f 
order —100 MeV or less attractive, as summarized recently in Ref. ?]]. For a recent overview 
of KN and ^-nucleus dynamics, see Ref. js]. 

The subject of multi-^ nuclei was studied in Refs. {3, 4], where the focal question con- 
sidered was whether or not kaon condensation could occur in strong-interaction self-bound 

n 

nuclear matter. Yamazaki et al. |9J argued that K mesons might provide the relevant phys- 
ical degrees of freedom for reaching high-density self-bound strange matter that could then 
be realized as multi-.K" nuclear matter. This scenario requires that beyond some thresh- 
old value of strangeness exceeds m^c 2 + /ijv — rn\c 2 > 320 MeV, where /ijv is the nucleon 
chemical potential, thus allowing for the conversion A — > K + N in matter. For this strong 
K binding, A and 5 hyperons would no longer combine with nucleons to compose the more 
conventional kaon-free form of strange hadronic matter, which is made out of {N, A, S} 
particle-stable configurations 10, H] (see Ref. [3] for an update), and K mesons would 
condense then macroscopically. However, our detailed calculations in Ref. J4] demonstrated 
a robust pattern of saturation for B^ and for nuclear densities upon increasing the number 
of K mesons embedded in the nuclear medium. For a wide range of phenomenologically al- 
lowed values of meson-field coupling constants compatible with assuming a deep ^-nucleus 
potential, the saturation values of B^ were found generally to be below 200 MeV, consider- 
ably short of the threshold value of ~ 320 MeV required for the onset of kaon condensation 

n 

under laboratory conditions. Similar results were subsequently published by Muto et al. [13| . 
Our discussion here concerns kaon condensation in self-bound systems, constrained by the 
strong interactions. It differs from discussions of kaon condensation in neutron stars where 
weak-interaction constraints are operative for any given value of density. For very recent 
works on kaon condensation in neutron-star matter, see Ref. jl^], where hyperon degrees 
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of freedom were disregarded, and Ref. 15(, where the interplay between kaon condensation 
and hyperons was studied, and references to earlier relevant work cited therein. 

In our calculations of multi-^ nuclei [4], the saturation of Br emerged for any boson-field 
composition that included the dominant vector w-meson field, using the F-type SU(3) value 
9ujKK ~ 3 associated with the leading-order Tomozawa- Weinberg term of the meson-baryon 
effective Lagrangian. This value is smaller than in any of the other commonly used models 
{4]. Moreover, the contribution of each one of the vector 0-meson and p- meson fields was 
found to be substantially repulsive for systems with a large number of antikaons, reducing 
Br as well as lowering the threshold value of the number of antikaons required for saturation 
to occur. We also verified that the saturation behavior of Br is qualitatively independent 
of the RMF model applied to the nucleonic sector. The onset of saturation was found to 
depend on the atomic number. Generally, the heavier the nucleus is, the more antikaons it 
takes to saturate their separation energies. We concluded that K mesons do not provide a 
constituent degree of freedom for self-bound strange dense matter. 

In the present work we extend our previous RMF calculations of muXti-K nuclei into 
the domain of multi--^ hypernuclei, to check whether a joint consideration of K mesons 
together with hyperons could bring new features or change our previous conclusions. This 
is the first RMF calculation that considers both K mesons and hyperons together within 



finite self-bound hadronic configurations. The effect of hyperonic strangeness in bul 



i on the 



dispersion of kaons and antikaons was considered by Schaffner and Mishustin 16) . More 



recently, kaon-condensed hypernuclei as highly dense self-bound objects have been studied 
by Muto [171 ] . using liquid-drop estimates. 

The plan of the article is as follows. In Sec. [Til we briefly outline the RMF methodology 
for multi-.K" hypernuclei and discuss the hyperon and K couplings to the meson fields used 
in the present work. Results of these RMF calculations for multi-A' hypernuclei are shown 
and discussed in Sec. IIHI We conclude with a brief summary and outlook in Sec. IIVI 
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II. MODEL 



A. RMF formalism 

In the present work, our interest is primarily aimed at multiply strange baryonic systems 
containing (anti)kaons. We employed the relativistic mean-field approach where the strong 
interactions among pointlike hadrons are mediated by effective mesonic degrees of freedom. 
In the following calculations we started from the Lagrangian density 

C = B [YfDp - {M B - g aB a - g a * B a*)\ B 

+ (D P K) ] (D"K) - {m\ - g aK m K a - g a * K m K a*)K^ K (1) 

+ (cr, cr*,u p , p p , M , A p free-field terms) - U(a) - V(u), 

which includes, in addition to the common isoscalar scalar (<r), isoscalar vector (u>), isovector 
vector (p), electromagnetic (A) fields, and nonlinear self-couplings U(a) and V(u), also 
hidden strangeness isoscalar a* and fields that couple exclusively to strangeness degrees 
of freedom. Vector fields are coupled to baryons B (nucleons, hyperons) and K mesons via 
the covariant derivative 

Dp = df, + i g^ uj p + i g^ I-p^ + i g^ <p^ + i e ( J 3 + \Y)A lt , (2) 

where <£> = B and K, with I denoting the isospin operator, J 3 being its z component, and 
Y standing for hypercharge. This particular choice of the coupling scheme for K~ mesons 
ensures the existence of a conserved Noether current, the timelike component of which can 
then be normalized to the number of K~ mesons in the medium, 

p K - = 2(E K - + g^Kid + g pK p + g (j)K (f) + eA)K + K', Jd 3 xp K -=K, (3) 

and serves as a dynamical source in the equations of motion for the boson fields in matter: 

(-V 2 + ml)a = g aB BB + g aK m K K + K- - ^-U(a) 

(-V 2 + mf)a*= g a * B BB + g a * K m K K + K~ 

d 

(-V 2 + ml )u = g wB B*B - g uK p K - + g^ v H (4) 
( -V 2 + m 2 )p = g pB BU z B - g pK pK- 
(- V 2 + ml)(f) = g^sB^B - g^KpK- 

-V 2 A = eB^(h + \Y)B-ep K -. 



These dynamically generated intermediate fields then enter the Dirac equation for baryons, 

[-ia • V + (3 (M B - g aB a - g a * B °*) + 9wB^ + 9 P Bhp + 9+b$ + e (h + \Y) A] B = eB 

(5) 

and the Klein-Gordon equation for K~ mesons, 

[- V 2 - E\- +m 2 K + U K -]K- = 0, (6) 
with the in-medium K~ self-energy, 

n^- = - gaKni K a - g a *Km K a* - 2E K - (g wK u + 9 p kP + 9^k<\> + eA) 
- {SukU + 9 p kP + 94>k4> + eA) . 

Hence, the presence of the K mesons modifies the scalar and vector mean fields entering the 
Dirac equation, consequently leading to a dynamical rearrangement of the baryon configu- 
rations and densities that, in turn, modify the K quasibound states in the medium. This 
requires a self-consistent solution of these coupled wave equations, a procedure followed nu- 
merically in the present as well as in our previous works. In the present work, for the sake 
of simplicity, we have suppressed the imaginary part of IIr-- arising from in-medium K~ 
absorption processes except for demonstrating its effect in one example. Note that, for the 
range of values B K - > 100 MeV mostly considered here, the effect of Im 11^- was found to 
be negligible (see Fig. 1 of Ref. [4j). 

B. Choice of the model parameters 

To parametrize the nucleonic part of the Lagrangian density ([1]) we considered the stan- 
dard RMF parameter sets NL-SH Q] and NL-TM1(2) Q, which have been successfully 
used in numerous calculations of various nuclear systems. 

In the case of hyperons the coupling constants to the vector fields were fixed using SU(6) 
symmetry. For A hyperons this leads to 

2 -y/2 

9uK = -g^N, g P A = 0, g<f,A = - g U N- (8) 

The coupling to the scalar a field, g^/ g aN = 0.6184 (0.623) for the NL-SH (NL-TM) RMF 



model, was then estimated by fitting to measured A-hypernuclear binding energies 20j. This 



essentially ensures the well depth of 28 MeV for A in nuclear matter. The coupling of the A 



TABLE I: K and K separation energies, B^ and B K - , respectively, calculated statically (in MeV) 
for a single antikaon Is state in several nuclei, using the NL-TM nuclear RMF parametrizations 
(TM2 for 12 C and 16 0, TM1 for 40 Ca and above) and vector SU(3) coupling constants, Eq. ([TO]) . 
The difference B K - — B^ is due to the K~ finite-size Coulomb potential. 





12 C 


16Q 


40 Ca 


90 Zr 


208pk 


Br 


44.8 


42.7 


49.8 


54.5 


53.6 


B k - 


49.0 


47.6 


59.2 


69.4 


76.6 



hyperon to the scalar a* field was fixed by fitting to the measured value ABaa ~ 1 MeV of 
the uniquely identified hypernucleus ^He [2l|. For S hyperons, SU(6) symmetry gives 

1 V2 
9u;S — ~9uiN, g P s — —g P N, SV>H — — 2— g^N- (9) 

Because there are no experimental data for H(A)-H interactions, we set g$~ = g a *~, = 
to avoid parameters that might lead to unphysical consequences and that, in addition, are 
expected to play a minor role (in analogy to the small effect, of order 1 MeV for B^-, 
found upon putting g^A and g a *A to zero, and as is demonstrated below in Fig. 7 within a 
different context). The coupling to the scalar a field was then constrained to yield an optical 
potential Re V=- = —14 MeV in the center of 12 C [221 ] . This corresponds to g a s = 0.299g CT Ar 
for the NL-TM2 RMF model. 

Finally, for the antikaon couplings to the vector meson fields we adopted a purely F-type, 
vector SU(3) symmetry: 

2g^ K = 2g pK = V2g^ K = 9fm = 6.04, (10) 

where g pn is due to the p — > 2tt decay width [7]. (Here we denoted by gyp the VPP 
electric coupling constant gypp-) Using this "minimal" set of coupling constants to establish 
correspondence with chirally based approaches, we calculate the single antikaon Is separation 
energies Bp^ and B K - listed in Table [H These separation energies are lower roughly by 
25 MeV than those anticipated from KN — Y>ix coupled-channel chiral approaches [7J , most 
likely because the K* vector-meson off-diagonal coupling is not included in the standard 
RMF formulation. The missing attraction, and beyond it, is incorporated here by coupling 
the antikaon to scalar fields a and a*. SU(3) symmetry is not of much help when fixing the 



coupling constants of scalar fields. Because there still is no consensus about the microscopic 



24J ] , in this work 



origin of the scalar a field and the strength of its coupling to K mesons 
we fitted g a x to several assumed K~ separation energies B K - in the range of 100 — 150 MeV 
for a single K~ meson in selected nuclei across the periodic table, as implied by the deep 
A^-nucleus potential phenomenology of Refs. . Furthermore, for use in multistrange 
configurations, the coupling constant to the a* field is taken from /o(980) — > KK decay to 
be g a *K = 2.65 [16j. The effect of the a* field was found generally to be minor. For a more 



comprehensive discussion of the issue of scalar couplings, see our previous work 



C. Inclusion of the SU(3) baryon octet 



We considered many-body systems consisting of the SU(3) octet N, A, E, and H baryons 



that can be made particle stable against strong interactions [10t 111! ]. The energy release 
Q values for various conversion reactions of the type BiB 2 — > B 3 B^ together with phe- 
nomenological guidance on hyperon-nucleus interactions suggest that only the conversions 
— > AA and S°n — > AA (for which Q ~ 20 MeV) can be overcome by binding effects. It 
becomes possible then to form particle-stable multi-{A^, A, 5} configurations for which the 
conversion EN — > AA is Pauli blocked owing to the A orbitals being filled up to the Fermi 
level. For composite configurations with E hyperons the energy release in the EA" — > AA^ 
conversion is too high (Q > 75 MeV) and, hence, it is unlikely for hypernuclear systems 
with E hyperons to be particle stable. 

III. RESULTS AND DISCUSSION 

In Refs. {3], Q] we studied multi-A' nuclei, observing that the calculated K~ separation 
energies as well as the nuclear densities saturate upon increasing the number of K~ mesons 
embedded dynamically in the nuclear medium. This saturation phenomenon, which is qual- 
itatively independent of the applied RMF model, emerged for any boson-field composition 
containing the dominant vector tu-meson field which acts repulsively between K mesons. Be- 
cause the calculated K~ separation energies did not exceed 200 MeV, for coupling-constant 
combinations designed to bind a single K~ meson in the range Bk- ~ 100 — 150 MeV, it was 
argued that kaon condensation is unlikely to occur in strong-interaction self-bound hadronic 
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matter. In this section we demonstrate that these conclusions hold also when adding, within 
particle-stable multistrange configurations, large numbers of hyperons to nuclei across the 
periodic table. 

A. Multi-{iV, A, K~} configurations 

Figure [H presents Is K~ separation energies Bk- in 16 + r/A + kK~ multi-fT'A hy- 
pernuclei as a function of the number k of K~ mesons for r] — 0, 2, 4, 6, and 8 A hyperons, 
calculated in the NL-SH model for two values of g a K (daK = 0.233^^ and 0.391p CT jv) chosen 
to produce B K - = 100 and 150 MeV, respectively, for r] = 0, k = 1. In addition, the lower 
group of curves with B K - < 60 MeV corresponds to g a K = 0. The figure illustrates satu- 
ration of Bk- with the number of antikaons in multi-A hypernuclei. There is an apparent 
increase of Bk- (up to 15%) when the first two A hyperons fill the Is shell. Further A 
hyperons, placed in the p shell, cause only insignificant variation of Bk- for small values of 
k. However, the effect of the lp3/2-shell hyperons increases with the number of antikaons, 
and for k = 8 it adds another 5 — 10 MeV to B^-- The separation energy B%- remains 
almost unaffected (or even decreases) by the next two A hyperons placed in the lpi/2 shell. 
The figure thus suggests saturation of the K~ separation energy also with the number r] 
of A hyperons in the nuclear medium. When the K~ coupling to the a field is switched 
off, g a K = 0, the K~ separation energy assumes relatively low values, B K - < 50 MeV, and 
decreases as a function of n when Im 11^- is considered (solid lines). In this case, the effect 
of K~ absorption is not negligible as illustrated by the dot-dashed line showing B^- for 
Im 11^- = 0. The effect of Im 11^-- ^ for B K - > 100 MeV in the upper groups of curves 
is negligible and is not shown here or in all subsequent figures. 

It is worth noting that 77 = 8 is the maximum number of A hyperons in our calculation that 
can be bound in the 16 nuclear core. In some of the 16 + rjA + kK~ allowed configurtions, 
lpi/2 neutrons became less bound than l<i 5 / 2 neutrons because of the strong spin-orbit 
interaction. (This occurs, e.g., for 77 = when k > 5 or for rj = 8 when k > 3.) However, 
the total binding energy of the system was found always to be higher for configurations with 
lpi/2 neutrons. Consequently, the standard shell configurations of oxygen are more bound 
and are thus energetically favorable. 

The saturation of B^- upon increasing the number of A hyperons in multi-i^'A hyper- 
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FIG. 1: (Color online) The Is K separation energy B K - in 16 O + r/A + kK as a function of 
the number n of antikaons for several values of the number rj of A hyperons, with initial values 
B K - = 100 and 150 MeV for 77 = 0, « = 1, calculated in the NL-SH RMF model. The solid 
(dot-dashed) lines with open symbols correspond to g a K = including (excluding) Im Hr- . 



nuclei based on a 16 nuclear core holds also when going over to heavier core nuclei. Figure [2] 
shows the Is K~ separation energy B K - in 208 Pb + rjA + kK~~ mvlti-K~A hypernuclei as 
a function of both the number k of K~ mesons and 77 of A hyperons, calculated in the 
NL-TM1 model for g aK = 0.133^ such that B K - = 100 MeV for 77 = 0, k = I. For any 
given number rj of A hyperons, Br- saturates with the number k of K~ mesons, reaching 
its maximum value for k — 12. Morever, B^- increases with the number of hyperons up to 
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FIG. 2: (Color online) The Is K separation energy B K - in 208 Pb + rjA + kK as a function 
of the number k of antikaons for several values of the number rj of A hyperons, with initial value 
B K - = 100 MeV for ry = 0, k = 1, calculated in the NL-TM1 RMF model. 

7] = 20, when it reaches its maximum value Bk- ^ 110 MeV for k = 12, and then starts to 
decrease with rj. Consequently, in the Pb configurations with 100 A hyperons and more than 
5 K~ mesons, K~ mesons are even less bound than in configurations with no A hyperons. 
The decrease of B K - with rj beyond rj = 20 is apparently related to a depletion of the central 
nuclear density in the presence of a massive number of hyperons in outer shells, as confirmed 
by some of the subsequent figures, because B K ~ is greatly affected by the central nuclear 
density. 
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B. Multi-{iV, A,E,K~} configurations 

When building up baryonic multi-{iV, A, S} configurations with maximum strangeness 
for selected core nuclei, we first started by filling up A hyperon single-particle states in a 
given nuclear core up to the A Fermi level. Subsequently, we added S° and S~ hyperons 
as long as the reaction [AN, r]A, fiE] — > [(A — l)N,rjA, (// — 1)2] + 2A was energetically 
forbidden (here, [...] denotes a bound configuration). Finally, we checked that the inverse 
reaction [AN, 77 A, fiEL] — > [(A + 1)N, (77 — 2)A, (// + 1)S] is kinematically blocked as well. 
These conditions guarantee that such {N, A, S} multistrange configurations are particle 
stable against strong interactions, decaying only via weak interactions. 

Clearly, the amount of 5 hyperons bound in a given system depends on the depth — V= of 
the S-nucleus potential. We adopted a value for g a ^ that gives V~ ira ~ c = Vjg+Vv = —18 MeV, 
corresponding to a depth of — Kf chr ' ~ 14 MeV for use in the Schroedinger equation 22]. 
For comparison, in some cases we also considered KP irac = —25 MeV. 

The 16 core can accommodate up to 77 = 8 A hyperons in particle-stable configurations, 
and the 16 + 8A system admits many more, of order 40 K~ mesons. However, we have 
not found any energetically favorable conversion A A — > "EN in 16 O + 77A + kK~ systems. 
Therefore, S hyperons are not part of any particle-stable multistrange configurations built 
upon the 16 core. While checking the energy balance in heavier systems with 40 Ca, 90 Zr, 
and 208 Pb nuclear cores, we found particle-stable configurations: 40 Ca + 20A + 2S°, 90 Zr + 
40A + 22° + 2H~, and 208 Pb + 106A + 82° + 18H~. We then embedded several K~ mesons in 
these configurations and studied density distributions and binding energies in such multi-K - 
hypernuclear systems. Figure [3] demonstrates the calculated Is K~ separation energy B K - in 
40 Ca+20A+2H°+/t J Fr-, 90 Zr+40A+2S°+2S-+fi;ir-, and 208 Pb+106A+8S°+18H~+fi;ir- as 
a function of the number k of K~ mesons. For comparison, in the case of the 208 Pb core, we 
also present calculations done excluding H hyperons but keeping the same number, r\ = 106, 
of A hyperons. A decrease of Br- upon adding hyperons (H in this case) is noted, in line 
with the trend observed and discussed for Fig. [2] above. 

The calculations shown in Fig. [3]were performed within the NL-TM1 nuclear RMF scheme 
using values of g aK = 0.211^ ( 40 Ca) and 0.163^ ( 90 Zr), which yield B K - = 100 MeV 
for a single K~ nuclear configuration with rj — /i = 0, where \x denotes the number of S 
hyperons. The figure demonstrates that the saturation of K~ separation energies, observed 
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FIG. 3: (Color online) The Is K separation energy B K - in 40 Ca, 90 Zr, and 208 Pb with t/A + ^H + 
kK~ as a function of the number k of antikaons, with initial value B K - = 100 MeV for rj = fj, = 0, 
k = 1, calculated in the NL-TM1 RMF model. 

for multi-A hypernuclei in Figs. [1] and [2], holds also when 5 hyperons are added dynamically 
within particle-stable configurations and that the heavier the system is, the larger number k 
of antikaons it takes to saturate B^- ■ It is worth noting that in all cases B^- does not exceed 
120 MeV. Finally, the two curves for a 90 Zr nuclear core in Fig. [3] (using diamond symbols) 
show the sensitivity to the value assumed for the H hyperon potential depth, the standard 
_y^)irac _ -^g Mey^ anc [ a somewhat increased depth — KP irac = 25 MeV, illustrating the 
tiny effect it exercises on B K - that is noticeable only for k < 12. 

A deeper H potential supports binding of more H hyperons in a given multi-A hypernu- 
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FIG. 4: (Color online) The Is K separation energy B K - in 40 Ca and 90 Zr with ^A+juH + K-fT , for 

yDirac = _ 25 M{J y 

clS cL function of the number k of antikaons, with initial value B K - = 100 MeV 
for rj = \i = 0, k = 1, calculated in the NL-TM1 RMF model. 

cleus. For V^ ira,c = —18 MeV, only 2S° and 2S° + 2H~ hyperons were found to be bound 
in 40 Ca + 20A and 90 Zr + 40A, respectively. However, for K? irac = -25 MeV it is possi- 
ble to accommodate up to 82° + 2H~ hyperons in Ca + 20A and 8H° + 8H~ hyperons in 
90 Zr + 40A. Figure H] presents the Is K~ separation energy B K - in multi-f^ - hypernuclei 
40 Ca + 20A + 8H° + 2H~ + kK~ and 90 Zr + 40A + 8S° + 8H" + kK~ as a function of the 
number k of K~ mesons, calculated in the NL-TM1 model for V? irac = —25 MeV, using 
values for g a x such that B^- = 100 MeV in 40 Ca + 1K~ and in 90 Zr + 1K~ . The figure 
illustrates that the saturation of the K~ separation energy occurs also in baryonic systems 
with three species of hyperons, A, S°, and S~, reaching quite large fractions of strangeness 
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[I-S'l/S = 0.57(0.8) for a Ca(Zr) core]. We note that the separation energy B K ~ barely 
exceeds 120 MeV in these cases too. 

We also studied the rearangement of nuclear systems induced by embedding hyperons 
and K~ mesons. Figure presents the evolution of the density distributions in Zr after first 
adding 40A + 4S hyperons (top panel) and then 10 K~ mesons (bottom panel). The nucleon 
density pjy in 90 Zr is denoted by a dotted line. The relatively weakly bound hyperons with 
extended density distributions (dashed line, solid diamonds) attract nucleons, thus depleting 
the central nucleon density p^ (dashed line, circles). Adding extra 10 K~~ mesons to the 
hypernuclear system induces large rearrangement of the baryons. The K~ mesons, which 
pile up near the origin (solid line, squares), attract the surrounding nucleons and hyperons. 
Consequently, the densities p^ and py (solid lines, solid circles and diamonds, respectively) 
increase considerably in the central region. The resulting configuration 90 Zr + 40A + 2H° + 
22~ + 10iT~ is thus significantly compressed, with central baryon density ps exceeding the 
nuclear density in 90 Zr by a factor of roughly 3. 

For comparison we present in Fig. [5] also the A hyperon (p\, open diamonds) and nucleon 
(Pn, open circles) density distributions calculated in 90 Zr + 40A + kK~ for k = and 10 K~ 
mesons. The removal of the ls-state H hyperons from the primary baryonic configuration 
90 Zr + 40A + 2S° + 2H~ affects considerably the hyperon density distribution py in the 
central region of the nucleus, this effect being magnified by the presence of K~ mesons. In 
contrast, the nucleon density pn remains almost intact. For k — 10, 5 hyperons appear to 
repel nucleons from the center of the multi-{A^, Y, K} system, much like A hyperons do. 



C. Multi-{iV, A,E,K+} configurations 

The -fT + -nucleus potential is known to be repulsive, with Vk+ ~ 30 MeV at central nuclear 
density Schaffner and Mishustin [l(| suggested that the presence of hyperons could lead 
eventually to a decrease of the repulsion that K + mesons undergo in nuclear matter so that 
the K + potential might even become attractive. Here we studied the possibility of binding 
K + mesons in hypernuclear matter, neglecting for simplicity dynamical effects arising from 
coupling K + mesons to the hypernuclear system. The i^ + -nucleus potential was constructed 
simply by applying a G-parity transformation to the corresponding K~ potential, choosing 
g a K such that it produces B K - = 100 MeV in the given core nucleus. 
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FIG. 5: (Color online) Density distributions in 90 Zr + 40A + 2~° + 2H~ + kK~ , for k = (top 
panel) and « = 10 (bottom panel), with B^- = 100 MeV for rj = \i = 0, k = 1, calculated in 
the NL-TM1 RMF model. The dotted line corresponds to the nucleon density pj\[ in 90 Zr. The 
densities pa (open diamonds) and ^at (open circles) in 90 Zr + 40A + Ki ; C _ are shown for comparison. 
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r(fm) 

FIG. 6: (Color online)The K + static potential in 16 + r]K - up, calculated in the NL-SH RMF 
model. 

Figure [6] shows the radial dependence of the real part of the static K + potential in various 
hypernuclear systems connected with 16 O. The dotted line shows the repulsive K + potential 
in 16 O for comparison. The figure indeed shows that the repulsion decreases, from roughly 
30 MeV down to roughly 20 MeV with the number of A hyperons added to the nuclear 
core, but the K + potential remains always repulsive in 16 0+i]A systems. Searching for a 
K + bound state in hadronic systems we also calculated the K + potential in more exotic 
multistrange hypernuclei A Z + rjA — up, where several protons are removed from the nuclear 
core in an attempt to increase the \S\/B ratio and to reduce Coulomb repulsion. Figure [6] 
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FIG. 7: (Color online) The K + static potential in 90 Zr+r/A+/x(H° + H"), calculated in the NL-TM1 
RMF model. 

indicates that such removal of protons from 16 O has a sizable effect on the shape of the 
K + potential, which may result in a shallow attractive pocket. However, the attraction is 
insufficient to bind a K + meson in these hadronic systems. Our calculations confirmed that 
the above conclusion holds also in heavier hypernuclear configurations based on Ca, Zr, and 
Pb cores. 

In heavier nuclei, where it becomes possible to accommodate also H hyperons in addition 
to A hyperons, the K + repulsion may be further reduced. This is demonstrated in Fig. [7] 
for a 90 Zr nuclear core. However, this reduction is insufficient to reverse the repulsion 
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into attraction. The figure also shows that the hidden strangeness couplings (chosen to be 
9iB — 2#iA, i = a*, 0) have no effect whatsoever on the reduction accomplished by the 
presence of H hyperons. 

Finally, we searched for K + bound states in nuclei sustained by K~ mesons. The presence 
of deeply bound K~ mesons makes the K + potential immensely deep (more than 100 MeV 
in 16 O + 8K~). Hovever, because the K~ mesons are concentrated at the very center of the 
nucleus, the K + potential is of a rather short range of about 1 fm. As a result, we found 
only very weakly bound K + states (by 1 MeV) in multi-{iV, Y, A - } configurations. A more 
careful treatment of K + K~ dynamics near threshold is necessary before coming to further 
conclusions, but our conclusion is not at odds with recent studies of the J = l/2,J 7r = l/2 + 



AAA system 
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IV. SUMMARY AND CONCLUSIONS 

In this work, the RMF equations of motion for multi-A hypernuclei were formulated and 
solved for self-bound finite multistrange configurations. The choice of coupling constants 
of the constituents - nucleons, hyperons, and A mesons - to the vector and scalar meson 
fields was guided by a combination of accepted models and by phenomenology. The sensi- 
tivity to particular chosen values was studied. The results of the RMF calculations show a 
robust pattern of binding-energy saturation for A mesons as a function of their number k. 



Compared to our previous RMF results for multi-A nuclei J4], the added hyperons do not 
bring about any quantitative change in the B K - («) saturating curve. The main reason for 
saturation remains the repulsion induced by the vector meson fields, primarily u, between K 
mesons. The SU(3)y values adopted here for g vK , Eq. (JTUJ), provide the "minimal" strength 
for g^K °ut °f several other choices made in the literature, implying that the saturation of 
Bk-(k) persists also for other choices of coupling-constant sets, as discussed in Ref. |4(. The 
repulsion between A mesons was also the primary reason for saturation in multi-A nuclei, 
both in our previous work and in Ref. [ljj. 

The saturation of Bk- with typical values below 200 MeV, considerably short of what 
it takes to replace a A hyperon by a nucleon and a A meson, means that A mesons do 
not compete favorably and thus cannot replace hyperons as constituents of strange hadronic 
matter. In other words, K mesons do not condense in self-bound hadronic matter. The 
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baryon densities of multi-A" hypernuclei are between (2 — 3)po, where po is nuclear-matter 
density. This is somewhat above the values obtained without K mesons, but still within the 
density range where hadronic models are likely to be applicable. 

Our conclusion of no "kaon condensation" is specific to self-bound finite hadronic systems 
run under strong-interaction constraints. It is not directly related to the Kaplan- Nelson con- 
jecture of macroscopic kaon condensation [271 ]. nor to hadronic systems evolving subject to 
weak- interaction constraints, such as neutron stars. Yet, this conclusion has been challenged 
recently by Muto [l^] who uses the liquid-drop approach to claim that multi-A' hypernuclei 
(termed by him "kaon-condensed hypernuclei" ) may provide the ground-state configuration 
of finite strange hadronic systems at densities about 9po- Of course this high value of density 
for kaon-condensed hypernuclei is beyond the range of applicability of hadronic models, be- 
cause quark-gluon degrees of freedom must enter in this density range. His calculation also 
reveals an isomeric multistrange hypernuclear state, without K mesons, at density about 2p 
which is close to what we find here within a RMF bound-state calculation. The appearance 
of a high-density kaon-condensed hypernuclear bound state in Muto's calculation might be 
just an artifact of the applied liquid-drop methodology, which does not provide an accurate 
substitute for a more microscopically oriented bound-state calculation. 

The role of K~ strong decays in hadronic matter was played down in the present cal- 
culation of multi-A"~ hypernuclei because our aim, primarily, was to discuss and compare 
(real) binding energies of strange hadronic matter with and without K~ mesons. The width 
of deeply bound K~ nuclear configurations was explored by us in Refs. , concluding 

that residual widths of order IV- ~ 50 MeV due to K~NN — > AN, UN pionless conversion 
reactions are expected in the relevant range of binding energy Bk- ~ 100 — 200 MeV. This 
estimate should hold also in multi-A"~ hypernuclei where added conversion channels are 
allowed: K~ NY -> AY, YY , K~NA -> NE, and K~ AY -> EY. We know of no physical 
mechanism capable of reducing substantially these widths, and therefore we do not antici- 
pate multi-A"~ nuclei or multi-A"~ hypernuclei to exist as relatively long-lived isomeric states 
of strange hadronic matter which consists of multi-{iV, A, E} configurations. 
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